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The elastic defoimation of two ellipsoidal solids in contact and subjected to a 
Hertzian stress distribution has been evaluated numerically as part of a general 
study of the elastic deformation of such solids in elastohydrodynamic contacts. 

In the analysis the contact zone is divided into equal rectangular areas and it is 
assumed that a uniform pressure is applied over each rectangular area. A study 
has been made of the influence on the size of the rectangular area upon accuracy. 
The results also indicate how far from the center of the contact one needs to go 
before elastic deformation becomes insignificant. 


SYMBOLS 

a = Semimajor axis of contact ellipse 
a = a/2m 

b = Semiminor axis of contact ellipse’ 
b = b/2m 

D = Defined by eq. (2) 

E == Modulus of elasticity 



F - Normal applied force 
h - Total film thickness 

h^ = Central film thickness due to elastohydro- 
dynamic lubrication 

m = Number of divisions of the semimajor or 
semiminir axis 

p = p/E\ dimensionless pressure 
p = Pressure 

^Professor, Dept, of Mechanical Engineering, 
Leeds Univ. , Leeds, England. 

P 

'Bearing Analyst, NASA Lewis Research Center, 
Cleveland, Ohio, U.S.A. 


R = Effective radius 
Rg = w/S 

r - Defined in fig. 3 

S = Film thickness due to the geometry of the 

solids, defined in eq. (8) 

W = F/E’R R , dimensionless load parameter 
X y 

w = Total elastic deformation 
w = Elastic deformation 

X^,X, X, xl Coordinate systems defined in the 
paper 

p = Poisson's ratio 
Subscripts: 

A = Refers to solid A 
B - Refers to solid B 

X, y = Refers to the coordinate system defined in 
the paper 

INTRODUCTION 

Elastohydrodynamic lubrication (ref. 1) is 
defined as the study of situations in which elastic 
deformation of the surrounding solids plays a sig- 



nificant role in the hydrodynamic lubrication proc- 
ess. This paper will not be concerned with the 
hydrodynamic lubrication process but only with 
the deformation due to pressure of one elastic 
solid upon another. Succeeding papers will deal 
with the cotqDling of the elastic deformation and 
hydrodynamic equations with the object of obtain- 
ing a complete three-dimensional pressure and 
film thickness profile. Reference 1 distinguishes 
between two forms of distortions which may exist 
in machine elements. The contact geometry may 
be affected by overall distortion of the elastic ma- 
chine element resulting from applied loads as 
shown in figure 1(a). 

In addition the normal stress distribution in 
the vicinity of the contact may produce local elas- 
tic deformations which are significant when com- 
pared with the lubricant film thickness, as shown 
in figure 1(b). This is the mode of deformation 
with which this report will concern itself. The 
important distinction is that the first form of de- 
formation is relatively insensitive to the distribu- 
tion and magnitude of the stresses in the contact 
zone, whereas the second mode of deformation is 
intimately linked to the local stress conditions. 

The deformation analysis will assume that the 
contact can be divided into rectangular areas and 
that the pressure can be assumed to be uniform 
within each rectangular area. Once the elastic 
deformation has been formulated, investigations 
will be performed to answer the following queries: 

(1) How fine do the semimajor and semiminor 
axes need to be divided to achieve a given accui> 
acy in deformation prediction? 

(2) How far out from the center of the contact 
is it necessary to go before deformation becomes 
insignificant compared with the natural separation 
of the solids ? 

These questions will be investigated for light and 
heavy applied loads, for equal spheres in contact, 
and for a contact that is common to the outer race 
of a ball bearing. The results of this investigation 
are given in greater detail in reference 2. 

GEOMETRY OF CONTACTING ELASTIC SOLIDS 

Two solids having different radii of curvature 
in a pair of principal planes (x and y) passing 


through the contact between the solids make con- 
tact at a single point under the condition of no ap- 
plied load. Such a condition is called point contact 
and is shown in figure 2 , Now when the two solids 
in figure 2 have a normal load applied to them, the 
result is that the point expands to an ellipse with 

being the semimajor axis and »^b^^ being the 
semiminor axis. The normal applied load in 
figure 2 lies along the axis which passes through 
the center of the solids and ihrou^ the point of the 
contact and is perpendicular to a plane which is 
tangential to both solids at the point of contact. 

The method used in evaluating the semimajor and 
semiminor axes is the same as that used in refer- 
ence 3 and therefore will not be repeated here. 

ELASTIC DEFORMATION 

Having the semimajor and semiminor axis of 
the contact ellipse the elastic deformation which 
occurs inside and outside the contact can be eval- 
uated. The approach to be used here will be one 
that can be quickly evaluated on the digital com- 
puter. The reason for this is that succeeding 
papers will deal with coupling the elastic deforma- 
tion with hydrodynamic equations thereby compli- 
cating things considerably. Therefore an elastic 
deformation analysis that is quickly evaluated and 
accurate will help assure success of succeeding 
papers of elastohydrodynamic lubrication of point 
contacts. Figure 3 shows a rectangular area of 
uniform pressure with the coordinate system to be 
used. From Timoshenko and Goodier (ref. 4) the 
elastic deformation at a point (X, Y) of a semi- 
infinite solid subjected to a pressure ^^p»^ at the 
point (X^, Yj^) can be written as 

2p dX^ dY^ 

dw = 

TrE’r 


The elastic deformation at a point (X, Y) due to the 
uniform pressure over the rectangular area 
2a X 2b is thus 
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where 


where 


P=J^ 

Integrating the above gives 

w = -PD 


where 


m = [k - i + 1 j 
n = |£ - j +l| 


(4) 

( 5 ) 


( 1 ) 


PRESSURE DISTRIBUTION 


D = (X+b)ln 
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/(Y+a)^+(X+b)^ 
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+ (Y + a)ln 


Within the contact ellipse the pressure will be 
assumed to be Hertzian. Therefore, using the 
coordinate system of figure 4 the dimensionless 
pressure is 


+ (X-b)ln 


(X+b) + 
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where 


W = - 




[(Y+a)+ V(Y+a)^+(X-b)^J 
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/ — 1 


N ^Ay ^By 

(X-b)+ ^ 

r(Y-a)^ + (X-b)2 

(2) 

The pressure outside the contact will be assumed 

(X+b)+ 1 

f (Y-a)^ + (X+b)2^ 


to be zero. 

in equation (1) represents the elas- 

FILM THICKNESS 


+ ( Y - a)ln 


tic deformation at a point (X, Y) due to a rectangu- 
lar area (2a x 2b) of uniform pressure ^’p. If 
the contact ellipse is divided into a number of 
equal rectangular areas, the total deformation at 
a point (X, Y) due to the contributions of the var- 
ious rectangular areas of imiform pressure in the 
contact ellipse can be evaluated numerically. 
Figure 4 shows a sample of dividing the area in- 
side and outside the contact into a number of equal 
rectangular areas. For purposes of illustration 
the contact was divided into a grid of 6 x 6 rec^- 
tangular areas. The effect of the fineness of this 
grid will be discussed later. Making use of fig- 
ure 4 the total elastic deformation at any point in- 
side or outside the contact ellipse due to the rec- 
tangular areas of uniform pressure within the con- 
tact can be written as 


The distance separating the two undistorted 
solids shown in figure 2 while using the coordinate 
system developed in figure 4 can be written as 


g - (X - b)^ ^ (Y - a)^ 
2R„ 


( 8 ) 


The total film thickness when a contact is elas- 
tohydrodynamically lubricated can be written as 

h = h^ +S(X,Y) +w(X,Y) (9) 

where 


Wi 


Kn 


> • • • 




i,j m,n ' ' 


central film thickness due to elastohydrody- 
namic lubrication 

w elastic deformation inside and outside the con- 
tact region 

The significance of the elastic deformation relative 
to the film thickness due to the geometry of the 
contacting solids can be e3q>ressed as 






INPUT CONDITIONS 

From figure 4 it can be seen that we need to 
concern ourselves with the following; 

(1) How fine must the divisions of and 

be? We will assume that the number of di- 
visions of and will be the same. There- 
fore, 

a b ^ V 

m = — = — (11) 

2a 2b 

In this paper we will let m = 3, 4, and 5. 

(2) How far from the semimajor and semi- 
minor axes must one go before R 2 (eq. (10)) be- 
comes insignificant? 

To check the accuracy of the elastic deforma- 
tion results for m = 3, 4, and 5 the number of 
equal divisions along the semimajor and semi- 
minor axes are increased by three times (m - 9, 
12, and 15) and then corresponding points are 
compared* The following equation describes the 
percentage accuracy of the results compared with 
the finest mesh size predictions. 



The limiting conditions that were evaluated on 
the computer are shown in table I. It was specu- 
lated that conclusions which could be made for 
these limiting conditions could also be made for 
any intermediate conditions. The four limiting 
conditions shown in table I are two extremes of 
applied normal load, a light load of 8.964 N (2 lb) 
and a heavy load of 896. 4 N (200 lb). The two 
extremes of curvature of the solids shown in 
table I are equal spheres in contact and a ball and 
outer race of a ball bearing. The elliptical eccen- 
tricity parameter (k - a/b) for the equal spheres 
in contact is equal to one and for the ball and outer 
ract it is equal to five. 


DISCUSSION OF RESULTS 

Tables n and III give the characteristics of the 
deformed shape of the contacting solids along the 
semimajor and semiminor axes when the axes are 


divided into five equal divisions and the dimension- 

-7‘ 

less load parameter is equated to 0,2102x10 and 
-5 

0.2102x10 . Some observations which can be 

made about these tables are: 

(1) Comparing table II with III, which amounts 
to changing the normal applied load from 8.964 N 
(2 lb) to 896.4 N (200 lb), the following can be con- 
cluded: 

(a) R 2 does not change in the correspond- 
ding tables. That is, regardless of the normal 
applied load, the ratio of the elastic deforma- 
tion to the natural separation of the solids is 
unchanged. 

(b) does not change in the correspond- 
ing tables. This condition is no doubt because 
of the condition mentioned in 1(a), 

(2) The separation due to the geometry of the 
contacting solids plus the elastic deformation 

(S + w) is very close to being constant in the con- 
tact, The value of (S + w) at the farthest point from 
the center of the contact and yet still in the contact 
differs the most from the other values of (S + w) in 
the contact, increasing slightly. 

(3) The percentabe difference in surface de- 
formation calculations for two mesh sizes differing 
hy a factor of three was shown to be small. For 
the worst case was found to be less than 2 per- 
cent. That is the elastic deformation for m == 5 at 
corresponding points and m = 15 differ by less 
than 2 percent, which is extremely good. 

(4) The ratio of the elastic deformation to the 
natural separation of the contacting solids, 

seen to decrease substantially as one moves away 
from the center of the contact. 

In order to better illustrate the results shown 
in the tables, figures 5 through 9 are presented. 

In figures 5 through 8 the solid curve represents 
the case of equal spheres in contact which is repre- 
sented by = By = 0. 558 centimeter (0. 2188 in. ) 
and the dotted curve represents the ball and outer 
race in contact which is represented by = 1.284 
centimeters (0.5055 in.), Ry = 15.00 centimeters 
(5,906 in.). Also due to the observation made in 
discussing the tables that R 2 and are not a 
function of the normal applied load the results 
shown in figures 5 through 8 apply for any normal 
applied load. 



Figure 5 shows the effects of the location 
along the semtmajor axis on the percentage differ-- 
ence in elastic deformation when m - 3 and 
m = 9. Here an *'edge effect^’ can be seen which 
is a rapid rise in percentage difference in the 
elastic deformation when m = 3 and for corre- 
sponding points when m = 9- This rapid rise is 
due to the pressure being either zero if the center 
of the rectangular area shown in figure 4 is out- 
side the contact or of order 10^ if the center of 
the rectangular area is within the contact* How- 
ever, it is speculated that in lubricated contacts 
where the pressure gradients are, in general, 
more gradual than those encountered near the edge 
of a dry Hertzian contact, this ^^edge effect^ ♦ is 
likely to be less significant. It is also to be noted 
in figure 5 that outside the contact, the value of 
Rg decreases rapidly. 

Figure 6 shows the effect of the location along 
the semimajor axis on the percentage difference 
in elastic deformation when m - 3, 4, and 5 and 
the more exact elastic deformation when m = 9, 

12, and 15, respectively. In this figure we see a 
large drop in R^ when going from m = 4 to 
m = 5. This also brings down the edge effect con- 
siderably, This therefore, makes a good case for 
letting m = 5 in any further computer evaluations. 

Figures 7 and 8 show the effect on the location 
along the semimajor and semiminor axes on the 
ratio of the elastic deformation to the distance 
separating the two solids in contact due to the 
geometry of the solids. Here we see the differ- 
ence that the shape of the contact has on how far 
from the semimajor and semiminor axes one needs 
to go before the elastic deformation becomes insig- 
nificant, To be more specific, from the curves we 
see that for equal spheres in contact (represented 
by solid lines in the figures) R 2 < 0. 05 corre- 
sponds to X > 2.6 b and y > 2.6 a. This means 
that the elastic deformation is less than 5 percent 
of the film thidkness due to the geometry effects 
when one moves away from the center of the con- 
tact a distance no less than 2.6 times the semi- 
major or Semiminor axes. For the ball and outer 
race in contact R 2 < 0,05 correspond to 
y>1.9a and x > 4, 0 b. This means that the 
elastic deformation is less than 5 percent of the 


film shape due to geometry effects when one moves 
only 1.9 times the semimajor axis away from the 
center of the contact and 4. 0 times the semiminor 
axis from the center of the contact. 

Figure 9 shows the effect of the location along 
the semimajor axis on the separation due to the 
geometry of the contacting solids plus the elastic 
deformation (S + w) when the dimensionless load 
parameter is 0. 2107x10”^, 0.5105x10 
0.2107x10'^, and 0.5105x10^^, This figure shows 
that (S + w) is constant within the contact. 

CONCLUSIONS 

A numerical analysis of the surface deforma- 
tion of two contacting ellipsoidal solids has been 
performed- The analysis assxmied that the pres- 
sure in the contact was Hertzian. It was assumed 
that the contact could be divided into rectangular 
areas with the pressure assumed to be uniform 
within each rectangular area. The resulting equa- 
tions were programmed on a digital computer. 

Four limiting conditions were evaluated on the 
computer. They consist of two extremes of applied 
normal loads, a li^t load of 8. 964 N (2 lb) and a 
heavy load of 896. 4 N (200 lb). The two other ex- 
tremes are of the curvature of the contacting 
solids. One of them is two equal spheres in con- 
tact and the other is that of a ball and outer race of 
a ball bearing. It was speculated that conclusions 
which could be made for the limiting conditions 
could also be made for any intermediate conditions. 

The results indicate that division of the semi- 
major and semiminor axes into five equal subdivi- 
sions is adequate to obtain accurate elastic defor- 
mation results. It was also found that the elastic 
deformation becomes insignificant compared with 
the normal surface separation for two equal 
spheres in contact at a distance from the center of 
2.6 times the semimajor axis. For a ball and 
outer race in contact it was found that a similar 
observation applied at a distance from the center 
of 1. 9 times the semimajor axis and 4. 0 times the 
semiminor axis. Finally, it was found that the 
separation due to the geometry of the contacting 
solids plus the elastic deformation (S + w) was al- 
most constant in the contact region. However, 
numerical values of (S + w) at points near the edge 



of the Hertzian contact show that a slight ’*edge 
effect" or error may be enco\jmtered in such re- 
gions, In lubricated contacts where the pressure 
gradients are, in general, more gradual than 
those encountered near the edge of a dty Hertzian 
contact, this effect is likely to be less significant. 
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(b) Local elastic deformation. 
Figure 1. - Types of elastic rieformatian. 
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Figure 2. - Geometry of contacting elastic solids. 



Coordinate, y 

Figure 5, - Effed of focatton a/ong the sernlmajor axis on the 
percentage difference in elastic deformation when m • 3 
and when m ” 9, 
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Figure 3. - Surface deformation of a semi -infinite body 
subjected to a uniform pressure over a reciangufar 
area. 
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figure 4. - Sample of how the area in and around the contact 
may be divided into equal rectangular areas. 
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Figure 6. ' Effect of locatfon elor»g me semimajor axis on the Figure 7. - Effect of the location along the semimajor axis on 

percentage difference in the elastic deformation when the ratio of the elastic detormatian to the distance 

m • 3, 4 and 5 and the more exact film shape when separating the two solids in contact due to the geometry of 

m • 9, 12 and 15, respectively, solids for m ^ 5. 
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Figure 8, - Effect of the location along the semiminor axis on 
the ratio of the elastic deformation to the distance separating 
the two solids in contact due to the geometry of the solids 
for m - 5, 
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Figure 9. - Effect of the location along the semimajor axis on 
the separation due to the geometry of the contacting solids 
plus the elastic deformation for m * 5. 



